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Abstract 

We present a meshfree method for the curvature estimation of membrane 
net works based on the Loca l Maximum Entropy approach recently presented 



in ( lArroyo and Ortid . l2006l ). A continuum regularization of the network is 
carried out by balancing the maximization of the information entropy cor- 
responding to the nodal data, with the minimization of the total width of 
the shape functions. The accuracy and convergence properties of the given 
curvature prediction procedure are assessed through numerical applications 
to benchmark problems, which include coarse grained mole cular dynamics 



simu l ations of the fluct uations of red blood cell membranes ( iMarcelli et al. 



20051 : iHale et al.l . l2009l ). We also provide an energetic discrete-to-continuum 
approach to the prediction of the zero-temperature bending rigidity of mem- 
brane networks, which is based on the integration of the local curvature 
estimates. The Local Maximum Entropy approach is easily applicable to the 
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continuum regularization of fluctuating membranes, and the prediction of 
membrane and bending elasticities of molecular dynamics models. 
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1. Introduction 



The estimation of the curvature tensor of membrane networks embed- 
ded in the 3D Euclidean space plays a key role in many relevant problems 
of diflFerential geometry, solid mechanics, biomechanics, biophysics and com- 
puter vision. Particularly important is the curvature estimation of fluctuat- 
ing bio-membranes, which are often modeled as particle networks, via molec- 
ular dynamics (MD) and/or coarse grained molecular dynamics (CGMD) 
approaches. The plasticity of cellular membranes is dependent on accu- 
rately selected mechanisms for sensing curvature and adopt different re- 
sponses according to the particular membrane curvature. These mecha- 
nisms depend on the interplay between proteins and lipids and can b e mod- 
ulated by cha nges in lipid composition (|Risselada and Marrinkl . l2009l ) , mem- 
brane fusion ( Martens and McMahon . 20081 ). formation of raft-like domains, 
oligomerization of scaflFolding proteins and/or insertion of wedge proteins 
into membranes. The dynamical chan ges in the membrane curvat ure can 
give rise to cell membrane remodelling ( jMcMahon and Gallop . 2005 ) result- 
ing in the formation of microenvironments that can facilitate the interac- 
tion of biomolecules in the cell. On a larger scale, these dynamical changes 
play a key role in controlling cellular growth, division and movement pro- 
cesses. Furthermore, as we already noticed, there has been a signiflcant 
amount of modeling work focussed on characterizing t he bending rigidity of 
ordered membranes (refer, e . g., to Seung and NelsonI ( |l988.. )). vesicle mem- 



branes (IGompper and Kro 



membrane ( IMarcelli et al 



1, 1996: 



Du et al 



2005; Hale et all 



2006h and the red blood cell 
20091 : bao 6^111 . l2006h . 



The continuum regularization of a membrane network is naturally per- 
formed through meshfree approximation schemes, which are well suited for 
the discrete-to-continuum scale bridging. Recently, a Local Maximum En- 
tropy (LME) approach has been proposed to construct smooth meshfree ap- 
proximants of given nodal data ( lArroyo and Ortid . l2006l : ICyron et al.l . l2009l : 
Li et al. 1, l20inh . The LME approach is a convex, non-interpolant approxi- 
mation method that suitably balances the maximization of the information 
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entropy corresponding to the given data (IJaines. Il957 ), with the minimiza- 
tion of the total width of the shape functions (IRajanl . 1 19941 ). Some of the 
distinctive features of such an approach consist of the non-negativity and the 
partition of unity properties of the shape functions, which in particular can 
be thought of as the elements of a discrete probability distribution; a weak in- 
terpolation (Kronecker-delta) property at the boundary ; and first- or h igher - 



order consistency conditions (I Arroyo and Ortiz . 2006 : Cyron et al. . 2009) 



. As compared to popular, 'explicit' approximation methods, such as ap- 
proaches utilizing B-Splines and Non-Uniform Rational B-Splines (NURBS), 
the LME requires more calculations and specifically the solution of a convex 
nonlinear optimization problem at each sampling point. Nevertheless, the 
LME guarantees hi ^h accuracy and smoothness of the continuous mapping 



(ICyron et al.l . |2009|), which are properties of fundamental importance when 
dealing with curvature estimation. Due to its mixed, local-global character, 
the LME approximation scheme can be conveniently used to filter the in- 
herent small scale roughness of the membrane, which is a distinctive feature 
of such an approach, as compared to popular computer graphic methods for 
estimating the curvatures of point-set surfaces, (e.g. moving least-squares 
(MLS) methods). Another peculiar advantage of the LME regularization 
consists of its ability to handle unstructured node sets, which do not re- 
quire any special pre-processing in such a scheme An extensive comparison 
of the application of the LME, ML S and B- Spline app roaches to structural 



vibration problems is presented in (ICyron et al.l . l2009l ) 



The present work deals with the formulation and the implementation 
of a curvature estimation method fo r membrane networks, w hich is based 
on the LME approach proposed in (lArroyo and Ortizl . l2006l ). To the au- 
thors' knowledge, such an application of the LME approximation has not 
appeared in the literature yet. In Section O we provide the explicit formulae 
for the second-order der ivatives of the LME shape functions (not given in 
flArrovo and Ortizl . l2006l )). and a LME procedure for the estimation of the 
lines of curvature and the principal curvatures at the generic node of a mem- 
brane network. Next, we present in Section [3] some numerical applications 
of the above procedure to membrane networks extracted from a sinusoidal 
surface and a spherical surface, establishing comparisons with exact solu- 
tions and assessing the convergence properties of the LME estimates. We 
also provide, in the same section, some numerical results abo ut the principal 



curvatures of the red blood cell (RBC) model proposed in (jMarcelli et al. 



20051 : iHale et al.l . l2009l ). and a discrete-to-continumm approach to the predic- 



3 



tion of the bending rigidity at zero temperature of MD membrane networks. 
The method and results presented here represent an essential first step to- 
wards an extensive estimation of the elastic moduli of the RBC, which will 
be the specific subject of future studies. Additionally, we plan to use the 
same approach to measure from trajectories of coarse-grain MD simulations 
the curvature of cell membranes affected by asymetric lipid bilayers or by 
protein-lipids interactions. Such future extensions of the present work are 
summarized in Section (U which also includes some final comments on the 
results presented in Section [31 



2. LME regularization of membrane networks 

2.1. Generalities on the LME approximation 

First, we will present how to find the continuum regularization of a 
given discrete set of nodes (or vertices) having Cartesian coordinates 
{xai, ^^2, Za} {a = 1, N ) with respect to a given frame {O, Xi, X2, z = X3}. 
We wish to construct a continuum surface described by the Monge chart 



N 



Z — Zn 



(x) = 5^ZaPa(x), 



(1) 



a=l 



where Pa are suitable shape functions of the position vector x = {xi,X2} in 
the x^ , x^, plane. When adopting the local maximum entropy approach pro- 
posed in ( lArroyo and Ortizl . I2OO6I ) , we determine the functions Pa by solving 
the following optimization problem: 



For a given x, we minimize: 

N N 
//3(X, p) = /3 |X - X^l^ + ^Pa logPa 



a=l 



a=l 



subject to: 



Pa>0, a = l,...,A^ 



(2) 



(3) 



N 



a=l 



(4) 
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N 

^p«x« = x (5) 

a=l 

where /3 G [0, +00) is a scalar parameter, and p = {pi, ...,PAr}. The con- 
straints (|3]), (|4j) and ([5j) enforce the non-negativity of the shape functions, 
the partition of unity property, and the first-order consistency conditions, re- 
spectively. It is worth noting that equations (|lj) and ([5j) g uarantee that affine 



funct i ons are exac t ly rep roduced by the LME scheme ( (I Arroyo and Ortid . 



eprc 

20061 : iCvron et all l2009h ). On the other hand, equations (|3j) and al- 
low us to regard p(x) as a discrete probability distribution, and the q uantity 
Hj{p ) = — J2a=i Pa log Pa as the corresponding information entropy ( (jjaines . 
19571 )). The quantity Vl^(p) = X^^iPa |x — x^^l^ instead represents the total 



width of the shape functions Pa at the given x. Depending on the value of 
the LME problem suitably balances the maximization of the information 
entropy corresponding to the given nodal data with the minimization of the 
t otal width of the shape functions Pa- A global maximum-entropy scheme 



(( iJainesl . I1957I )) is recove red by settin g /3 = in ([2]), while a minimum- width 
approximation scheme (( Raj an . 19941 )) is obtained in the limit /3 +oc. 

Now, we introduce the partition function Z(x, A) = ^^-^^^(x. A), 
where A = {Ai,A2} denotes the vector of the Lagrange multipliers of the 
first-order consistency conditions ([5j), and it results in 



Z«(x, A) =exp [-/3|x-x^| + A • (x - x«)] . (6) 

It can be shown (( j Arroyo and Ortiz . 20061 )) that, for any x G convX, the 
LME problem admits the unique solution shown below 



p:(x) 



^a(x, A*) 

Z(x,A) 



(7) 



where 



A* = arg min {F(A) = log Z(x, A)} (8) 

It is useful to the employ the Newton-Raphson method to solve equation 
dH]) iteratively. Let A^ denote the approximate solution to (|8]) at the kth 
iteration. A Newton-Raphson update furnishes 
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(9) 



where and 3^ are the particularization of the gr adient r and the Hessian J 
of F for A = A^. Straightforward calculations (cf. ( j Arroyo and Ortiz . 20061 )) 



give 



r(x, A) = VF(A) = J2 = T.P^(^ " 

a=l a=l 



(10) 



J(x,A) = V'F(A) 



2-^ fl\ 



a=l 
N 



(x - X^) 



= ^pa(x-x«)(8)(x-x«)-r®r (11) 



a=l 



2.2. Derivatives of the LME shape functions 



The analysis carried out in (lArroyo and Ortizl . l2006l ) leads to the following 
expression of the spatial gradient of 



VK = -K(J*)-^(x-x„) 



(12) 



where J* = J(x, A*). 

Now we compute the second-order derivatives of the LME shape func- 
tions. Differentiating both sides of (fT2l) once, we get 



P. 



dxidxj 



_d_ 



The only quantity that needs to be computed on the right-hand side of 
(fT3l) is Ji^j- Differentiating both sides of the identity shown below 



T* T*-l ^ 



(14) 



where 6j denotes the Kronecker symbol, leads to 



T* — 1 



T* — 1 T* T* — 1 

'^im mn,j nk ' 



On the other hand, from ( fTTI) we deduce the result 



N 



a=l 



~ N 

-1 

~ N 



a=l 



a=l 



(15) 



(16) 



Equations (fT5l) and (fT6l) allow us to derive the explicit formulae (fT3l) . 



^.5. Lme5 of curvature and principal curvatures of membrane networks 

Let us now examine a node set = {{xai^Xa^^ Za}^ a = 1, A^} ex- 
tracted from a membrane network lying in the 3D Euclidean space (Fig. [1]). 
The Monge chart 



N 



(17) 



defines the LME regularization of that we will denote by Sn in the 
following. The unit vectors 1/(1), 1/(2} are tangent to the lines of curvature of 
Sn^ and the principal curvatures /ci,A;2 of su ch a su r face c o rrespond to th e 
solution of the eigenvalue problem (see, e.g., (jStokerl . 1 19691 ). (INaghdil . 1 19721 ): 
Appendix A. 2) 



/^tM ^er) =0 (7=1,2) 



(18) 



where aa(3 and ha(3 are the first and the second fundamental forms of Sn^ 
defined by 
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By taking into account the previous expressions of the shape functions 
]9* and their derivatives, we easily derive from Eqns. f lTSl) and ( IT9l) the 
LME estimates of the mean curvature Hf^ = l/2(/ci + /C2), and the Gaussian 
curvature = ^1^2 of at the given x. 



Figure 1: (Color online) Node set extracted from a 3D membrane network and its orthog- 
onal projection onto the xi^X2 plane. 



3. Numerical results 

3.1. Sinusoidal membrane 

We begin by considering the node set such that and are 
randomly generated numbers within the interval [0, tt], and it results in Za = 
sin(x^^ + Xa2)- Fig. [2] shows the LME surfaces Sn obtained in this case 
for = 500, /3 = 0.001, and /3 = 10. The ^^v are sampled over a 12 x 12 
uniform grid of points defined over the Xi, X2 region [0, tt] x [0, tt]. We observe 
from Fig. [2] that the Sn corresponding to /3 = 0.001 is almost fiat, while the 
Sn corresponding to /3 = 10 fairly reproduces the local shape of Xn in the 
neighborhood of each node. 



8 



P = 0.001 



Figure 2: (Color online) LME approximations of scattered data extracted from the sinu- 
soidal surface S : z = sin(xi^ + X2) for two different values of /3. 

Next, we examine uniform grids of nodes = {^01,^02} ^^^^ the Xi,X2 
region D = [0,3] x [0,3], and the 3D node set = {x^ = {^^a^^a = 
sm{xl^ + Xa^)}, a = 1, A^}. For each x« G L^' = [0.5, 2.5] x [0.5, 2.5] C D, 
we further consider the subset C X^^ which is generated by the mth 
nearest neighbors of x^, m being an integer parameter. We employ the 
node set X^ to get LME estimates and of the mean and Gaussian 
curvatures of X^ at x^- Iii order to normalize the effects of /3 on the LME 
estimates, we rescale such a parameter as follows 



where /3 is a dimensionless quantity and it results in 

diam(X^) = max {|x«-X5|}. (21) 

It is useful to compare and with the 'exact' counterparts and 
K^^ which are easily computed through (fT8l) and (fT9l) . by replacing zn with 



(diam(X^))' 



(20) 
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z = sin(xi^+X2). The accuracy of the nodal LME estimates Hf^ = {H^^ a = 
1, A^'} and = {^at, ^ = 1, A^'} can be inspected by examining the 
foUowing Root Mean Square Deviations (RMSD) 



i (22) 

RMSDiKr,) = {^tAK% - K-f) /N' 

for different values of ^ and m. Here, A^' denotes the total number of nodes 
belonging to D^ We examine the following three different mesh sizes: h = 
0.0566 (A^ = 54 X 54); /i = 0.0405 (A^ = 75 x 75); and h = 0.0303 (A^ = 
100 X 100). Here and in the following examples, we solve the nonlinear 
optimization problem ([8j) using recursive Newton- Raphson updates ([9j), up 
to the termination condition |r^| < 10~^ diam(X^). Fig. [3] illustrates how 
the quantities RMSD[Hn) and RMSD{Kn) vary with ^ and h for fixed 
m = 10, while Figs. [H [5] and [6] depict 2D and 3D density plots of the data 
sets Hn and for several values of m, keeping ^ = 150, and h = 0.0303 
fixed. The results shown in Figures [H [5] and [6] point out that the LME 
estimates and exhibit uniform asymptotic convergence to the exact 
solutions H and respectively, for ^ > 100, m > 9, and h < 0.0303. 




Figure 3: (Color online) Root Mean Square Deviations (RMSD) of the LME approxima- 
tions to the mean curvature H and the Gaussian curvature K of the sinusoidal surface 
z = sin{xl + X2) over the xi, X2 domain [0.5, 2.5] x [0.5, 2.5]. Left: RMSD as a function of 
for m = 10 and = 0.0303. Right: RMSD as a function of h, for m = 10 and = 150. 
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3.2. Spherical membrane 

We examine, in the present and the foUowing sections, a closed membrane 
net work Xm correspond i ng to the CGMD model of the red blood cell given 



in ( iMarcelli et al.L l2005l : iHale et al.L l2009l ) . Such a model describes the sys- 



tem formed by the cytoskeleton spectrin network, the lipid bilayer and the 
transmembrane proteins of an actual RBC membrane, through a network of 
(virtual) particles embedded in a closed polyhedral surface showing M tri- 
angular facets. Each particle has sixfold coordination, with the exception 
of twelve 'defects,' which instead have fivefold coordination. The particles 
represent discrete areas of the RBC membrane and their equilibrium dis- 
tance ro is set equal to the average length of the spectrin filaments (^100 
nm). Each particle has mass m and is connected to it's nearest neighbors 
though linear springs of stiflFness fc, which are parameterized such that the 
network of particles reproduces the membrane rigidity of a RBC. The bend- 
ing rigidity is also accounted for by introducing dihedral angle potential 
energy terms of angular stiffness D between adjacent triangles. The global 
surface area of the polyhedral membrane is kept constant using Lagrange 
multipliers, in order to mimic the relative incompressibility of the lipid bi- 
layer of a normal RBC. The model under consideration is also able to keep 
the volume enveloped by the polyhedral membrane constant, with the aim 
to resemble the typical behaviour of a RBC in normal conditions. The 
dynamics of the RBC model can be obtained by integrating the Newton 
equations of motion o f each virtual partic l e usin R a standard MD code. We 
used DL POLY 2.20f flSmith and Foresteil . Il999l )). and employed the Nose- 
Hoover thermostat; 6 x 10^ steps with time-step At = 2.07 x 10"^ to, where 
to = ^m/k] N = 5762 particles; m = 5.82625 x IQ-^o kg; k = 8.3 /xN/m; 
D = 130 X 10~^° J; constant absolute temperature T = 309 K; constant sur- 
face area A = 4.986 x 10'' nm^; and constant volume V = 3.311 x 10^° nm^. 
Such parameter settings allow the LME re gularization to appro ximate a 



spherical surface S with radius r = 1992 nm ( IMarcelli et al.L 120051 ) 



The theoretical average surface S clearly has uniform principal curvatures 
ki = k2 = H = —1/f = —5.019 X 10~^ nm"-*^ (the minus sign follows from 
the outward orientation of the normal vector), and Gaussian curvature K = 
kik2 = 25.19 X 10"^ nm-2. Nodal LME estimates Hn = {H^, a = 1, A^} 
and Kn = {^at, cl = 1, A^} of the mean and Gaussian curvatures of the 

network can be obtained by introducing different local frames 5 ^i, ^2, ^} 
at each different node x^^, with Xi, X2 and z tangent to the local parallel, 
meridian, and radial lines, respectively. Table [1] shows the mean values and 
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the standard deviations of the data sets and defined as foUows 



Hn = llNY.a=iH%, sa{H^) = ^(e1i M - ^iv)^) /(iV - 1) 

(23) 



(24) 

for /3 = 100 and diflFerent values of m. The results in Table [1] highlight a 
good agreement between LSM estimates and exact solutions for m > 10. It 
has to be considered that the MD model doesn't reach a perfectly (average) 
spherical shape at equilibrium, due to the presence of the fivefold defects. 





m = 5 


m = 7 


m = 9 


m = 10 


m = 12 


Hn X 10^ a ' 


39.0626 


-7.1062 


-5.1293 


-5.0443 


-5.0311 


Sd{HN) X 10^ A"' 


2.8047 


0.1467 


0.0296 


0.0271 


0.0238 


Kn X IQio A"' 


1540.19 


50.488 


26.3077 


25.4453 


25.3118 


SdiKN) X IQio A"' 


234.04 


2.1551 


0.3152 


0.2840 


0.2473 



Table 1: Mean and standard deviation of the data sets Hn and Kn for a spherical mem- 
brane network with N = 5762 nodes and radius r = 1992 nm, considering P = 100 and 
different values of m. Exact solution: H = -5.019 x 10"^ A"\ K = 25.19 x 10"^° A"^. 



3.3. Principal curvatures of the RBC membrane 

We analyze in the present section a sl i ght different form ulation of the 
RBC model given in ( jMarcelli et al.l . l2005l : iHale et al.l . l2009l ). which diflFers 
from that discussed in the previous section only in terms of the volume 
constraint. Here, we set the volume enveloped by the RBC membrane equal 
to 0.65 times the volume of the sphere analyzed in the previous example, 
allowi ng the current mo del to assume the typical biconcave shape of a normal 
RBC ( IHale et al.l . 120091 ). The surface area of the RBC membrane is again set 
to 4.986 X 10'' nm^, as in the previous case. We compute nodal LME 



estimates H 



N 



a 



1,..., A^} and Kn = {K] 



1, 



A^} of 



the mean and Gaussian curvatures of a real RBC membrane by processing 
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the rolling average configurations of the CGMD model up to different 
simulation times t. Let n denote the weighted unit normal to the current 
vertex of the tri anRulation ass ociated with X^^ assuming the triangle 



areas as the weights (ITaubinl . 1 19951 ). In the present case, we define the local 
Xi axis as the direction of the edge attached to that has the minimum 
deviation from the parallel drawn on an ideal sphere passing through the 
same point. In addition, we define X2 by means of the vector product of the 
unit vectors in the directions of n and Xi, and z via the vector product of 
the unit vectors in the directions of Xi and X2. A graphical representation 
of the local parameterization introduced for the analyzed model of the RBC 
membrane is provided in Fig. [71 

Fig. [8] shows the dependence of the LME estimates and K^^^ = K^A 
on ^, for t = 75.5 to and m = 10. It is worth noting that represents an 
estimate of the total curvature K^^^ = K dA of the RBC model. Due to 
the 'Gauss-Bonet theorem' such a quantity only depend s on the Renus of Sn 



and should be equal to 47r, as in the case of a sphere ( Stoker . 1 19691 ). One 



observes from Fig. [8] slight oscillations of and Kj^^ with /3; Kj^^/A ^ 3.00 
for ^ < 100; and ^ 3.14 for ^ = 125. 3D density plots of the data sets 

Hjsf and are given in Figs. [91 and [TOl considering ^ = 125, m = 10, and 
different simulation times t. The results shown in Figs. [9] and [10] indicate 
that the LME regularizations of the average CGMD configurations smoothly 
describe the geometry of a normal RBC membrane. The LME regularization 
is indeed able to reproduce the biconcave shape of such a membrane (consider 
that the hidden bottom edges of the surfaces shown in Figs. M and [TOl are 
nearly specular with respect to the top ones), furnishing positive (red) and 
negative (blue) mean curvatures in correspondence with concave and convex 
portions, respectively (cf. Fig. [9]), and negative (blue) Gaussian curvatures 
over saddle-shaped regions (cf. Fig. [TOj). One observes perturbations in 
the mean and Gaussian curvature maps in correspondence with the fivefold 
defects of the membrane triangulation during the initial phase of the MD 
simulation. The MD simulation and the LME regularization are however 
able to smooth out the local noise produced by such perturbations, as the 
simulation time progressively increases. 

3.4' Estimation of the bending rigidity of membrane networks from MD sim- 
ulations 

Let us consider now a triangulated membrane network endowed with the 
following dihedral angle energy 
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E 



idihedral 



neighbors 



1 - cos A A' 



Fa - 



A,A^eCA 
neighbors 



(25) 



where Ca denotes the set of aU triangles forming the network; A A'' is the 
dihedral angle between the triangles A and A''; tt^a is the unit normal to the 
triangle A; and the summation runs over all the pairs A . A^ g Ca which sha re 



a com mon side. According to ( ISeung and Nelsonl . 1 19881 : iGompper and Krolll . 
19961 ) , the l imiting bending energy of su ch a network is a Helfrich-type bend- 
ing energy ( jPeuling and Helfrichl . 1 19761 ) of the form 



E 



bend 



^ I {{2Hf - 2K) dS 



(26) 



where k^h is the bending rigidity, and it results i^m = \/3D/ 3 for a sphere, 
and HiM = a/3D/2 for a c ylinder (Gompper and Kroh . 19961 ) (Note that in 
( Seung and Nelson . 1988 : Gompper and Krolj 1996l )the symbol H is used 
to denote twice the mean curvature of the network) . 

The LME prediction of the principal curvatures of the network allows 
us to estimate the limiting bending energy (l26l) . As a matter of fact, the 
computing of the mean curvatures = {H^^ a = 1, A^} of the current 
configuration leads us to approximate (l26l) as follows 



= — 



2 



N 



J2 - Svr 



(27) 



where is the surface area of the a-th element of a dual tessellation of 
the network, which we assume is formed by polygons joining the midpoints 
of the network edges with the triangle barycenters (barycentric dual mesh). 
Equation ( l27l ) can be employed to estimate the bending modulus at zero tem- 
perature of the network, t hrough the Cauchy-Born contribution ( Ericksenl . 
20081 : IZhou and Joosl . Il996l ) to the isothermal bending rigidity, here denoted 
by tv^. The latter consists of the configurational average of the instanta- 
neous bending rigidities of the fiuctuating network. On matching E^^^ to 
^dihedral (^eucrgctic discretc-to-continuum approach), we identify the instan- 
taneous bending rigidity with the ratio ^'^ihedray ^^bend^ L^^-i ' compute 
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/^^ through 

/ E^idihedral 



where (•) denotes the configurational average symbol. 

For the sake of example, let us reconsider now the MD model of the RBC 
membrane examined in the previous section. Fig. [TT] shows the time his- 
tory that we obtained for the zero-temperature bending rigidity /^^ of such a 



model, on computing the rolling averages of the quantity 2 ^'dihedray^bend 
during a MD simulation (cf. Sect. 13.31) . It is seen from Fig. [TT] that 
the value of /^^ slightly fluctuates during the simulation, featuring oscil- 
lations with progressively smaller amplitude as the computational time t 
increases. We estimated a limiting value of /^^ approximatively equal to 79 
J, for t > 100 to. Such a val ue is just slightly greater than that predicted 
in ( IGompper and Kroli 119961 ) for the sphere {^/3D/3 = 75.06 J), which is 
not surprising, since the biconcave (average) shape of the RBC model under 
consideration has the same genus of the sphere. 



4. Concluding remarks 



We have presented and numerically tested a meshfree approach to the 
curvature and bending rigidity estimation of membrane networks, through 
a sui table extension of the LME method formulated in ( lArroyo and Orti zl . 
20061 ) . The results of Sections 13.11 and 13.21 demonstrate the convergence 
properties of the LME curvature estimates, both for a rectangular geometry 
(flxed {xi, X2, z} axes), and for a spherical membrane. On the other hand, 
the results presented in Sections 13.31 and 13. 4[ emphasize the ability of the 
LME approach in tracking the local curvatures a nd the bending rigidity of the 



RBC model presented in ( jMarcelli et al.l . l2005l ) . Concerning the parameter 



estimation, we have found that the limitations P > 100 and m > 10 generally 
ensure stable curvature predictions. In particular, for m > 10 we found 
that 100 is approximately the smallest value of ^ that guarantees asymptotic 
convergence of LME curvature predictions for the examined examples. The 
value of ^ rules the degree of locality of LME approximations, which reduce 
to piecewise afline shape functions supported by a De launay triangulation 



for /3 ^ oc flArrovo and Ortizl . l2006l : ICvron et al.l . [20091) 



In closing, we suggest some directions for relevant extensions of the present 
work. We intend to apply the LME regularization algorithm to predict the 
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entire set of the isothermal membrane and bending rigidities of fluctuating 
biomembranes modeled through MD simulations, with special reference to 
the RBC membrane. W e plan to comb i ne th e LME regularization of the 



RBC model proposed in (IMarcelli et all 120051). with the elastic moduli es- 



timation procedures Riven i n (ISeung and Nelsonl . Il988l : iGompper and Kroli 



19961 : IZhou and Joosl . 119961 ). Another application of the LME procedure 



presented in this work is to mesure the curvature of lipid bilayers as mod- 
elled with MD simulations. Lipid bilayers are generally flat, however several 
diflFerent protein-driven processes will result in the membrane curvature that 
is required for various cell processes (i.e. fusion). One mechanism is caused by 
protein domains inserting shallowly into one of the lipid leaflets, which push 
the neighboring lipid head groups away and therefore causing local sponta- 
neous curvature. We are currently conducting coarse-grain MD simulations 
that model the interactions between diflFerent varieties of these protein do- 
mains and lipid bilayers. The LME procedure will allow us to quantify the 
amount of curvature that results from these interactions. 
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Figure 4: (Color online) 2D density plots of the LME approximations to the mean curva- 
ture H of the sinusoidal surface z = sin{xl-\-X2) over the xi, X2 domain [0.5, 2.5] x [0.5, 2.5], 
for ^ = 150, h = 0.0303, and different values of m (lower bound of the color bar: -4; .upper 
bound: +4). 
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Figure 5: (Color online) 2D density plots of the LME approximations to the Gaussian 
curvature K of the sinusoidal surface z = sin{x\ + X2) over the xi^X2 domain [0.5,2.5] x 
[0.5, 2.5], for P = 150, h = 0.0303, and different values of m (lower bound of the color bar: 
-0.56; upper bound: +0.56). 
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Figure 6: (Color online) 3D density plots of the LME approximations to the mean and 
the Gaussian curvatures of the sinusoidal surface z = sin{x\ + X2) over the X\^X2 domain 
[0.5, 2.5] X [0.5, 2.5], for ^ = 150, h = 0.0303, and m = 12. 
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Figure 7: (Color online) 3D map of the local bases {x^, xi, X2, z} introduced at selected 
nodes of a triangulated model of the RBC membrane {xi: yellow, X2' blue, z: cyan) . 
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Figure 9: (Color online) 3D density plots of the LME approximation s t o the mean 
curvature of of a CGMD model of the RBC membrane iMarcelli et al.l (|2QQ5[ ): lHale et al 
(|2QQ9h at different simulation times t, for ^ = 125 and m = 10 {Hn = Hn x 10^ A"^ 
lower bound of the color bar: — 11.5xlO~^A ; upper bound: +2.5 x IQ-^ A ). 
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Figure 10: (Color online) 3D density plots of the LME appro ximations Km to th e 
Gaussian curvatu re K of of a CGMD model of the RBC membrane iMarcelli et al.l (|2005[ ): 
Hale et al. (|2009[ ) at different simulation times for ^ = 125 and m = 10 (lower bound of 
the color bar: —2 x 10~^ A ; upper bound: +7 x lO'^ A ). 
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Figure 11: LME estiraates of the zero-teraperature bending rigidity of the RBC membrane 
model presented in Section [331 ^tt different times of a MD simulation. 
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